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Abstract 

This paper shows that there are no physical walls in the deconfined, high- 
temperature phase of Z(2) lattice gauge theory. In a Hamiltonian formulation, 
the interface in the Wilson lines is not physical. The line interface and its 
energy are interpreted in terms of physical variables. They are associated 
with a difference between two partition functions. One includes only the 
configurations with even flux across the interface. The other is restricted to 
odd flux. Also, with matter present, there is no physical metastable state. 
However, the free energy is lowered by the matter. This effect is described in 

terms of physical variables. 
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I. INTRODUCTION 



This paper discusses some aspects of the global Z(N) symmetry [l] of finite-temperature 
SU(N) or Z(N) gauge theory. It contributes to recent discussions of the physics of Z(N) 
phases and interfaces [0-0]. The example of Z(2) gauge theory is treated in detail here. 
Before focusing on that case, a more general introduction is given. 

There is a low-temperature phase of the theory, which is presumed to confine fundamental 
sources, and a high-temperature phase, which does not. Without matter fields, there is an 
order parameter — the finite-temperature expectation value of the Wilson line (L). The 
fundamental-representation line L(i) is the normalized trace of the product of the group 
elements on the links in the inverse-temperature direction at spatial site i. 

(L) = (L®) = (irTr\RU(i,U,4))) (1-1) 

The line carries a nontrivial representation of the global Z(N) symmetry. In the confining 
phase, (L) = 0, and the ensemble is Z(N) symmetric. In the high T phase, (L) takes 
one of N distinct values proportional to the Nth roots of unity z in Z(N), and the Z(N) 
symmetry is broken. The broken-symmetry, pure phases appear in Monte Carlo simulations, 
where runs of finite duration in the high-temperature phase and not too close to the critical 
temperature show the Wilson line fluctuating around one of the N available values. 

In the Hamiltonian description ||, the physical variables are the group elements on the 
links of the spatial lattice. In a Lagrangian formulation, there are also group elements 
on links in the inverse-temperature direction. These are unphysical, auxiliary variables 
introduced to enforce the Gauss law constraints. The Wilson line is constructed from the 
unphysical variables. It is a projection operator that forces the gauge field to be in a 
fundamental rather than a singlet state at the spatial position of the line. The global Z(N) 
symmetry of the Lagrangian formulation is not physical; it acts as the identity on all physical 
states 0-^,0]. There is a single physical, high-temperature phase, which is the same for all z. 
In the Hamiltonian formulation, the high-temperature phase is not distinguished by physical 
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broken symmetry. Rather the high-temperature phase has a percolating flux network that 
is not present in the low-temperature phase. 

Now consider Z(N) walls. If there were N, degenerate, physically-distinct, high- 
temperature phases, then there could be large regions of space in different phases and 
separated by interfaces. These walls could have a physically-significant interface tension. 
In the line variables, there are phase separations and interfaces. This topic has received 
much attention in the literature and the interface tension for the lines has been mea- 

sured in Monte Carlo simulations |§ . However, as noted above, the N phases of the Wilson 
lines are not physically distinct. Thus one may question the physical relevance of the line 
interfaces. There is work in the literature which is related this point and which asserts that 
the interfaces are not physical |4]|| . There is considerable overlap of the present paper with 
some of that work. Some of the statements in those papers are confirmed or derived by 
other means here. Other issues treated here have not been previously addressed. 

I use the Hamiltonian description of finite-temperature gauge theory on a lattice in three 
spatial dimensions in the A = gauge with Gauss law constraints. In this formulation, the 
physical interpretation is most accessible. I refer to this as the physical theory. The physical 
variables are the group elements on the links of the spatial lattice. They may be described 
with definite position on the group or in a basis of definite irreducible representation matrix 
elements, which is called the flux basis. Physical states are configurations of physical vari- 
ables subject to Gauss law constraints. With that definition, the term "physical" takes on 
a specific, technical meaning. Alternative formulations of the quantum gauge theory may 
lead to different pictures. 

The constraints are enforced through the use of additional group elements on sites. 
These are unphysical, auxiliary variables. In a four-dimensional, Lagrangian formulation, 
they become the elements on links in the inverse-temperature, 4-direction. The Wilson line 
is constructed from them. 

The interface exits in the line variables. It can be created and studied with the application 
of appropriate boundary conditions. Applying the boundary conditions and using physical 
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variables gives a description in terms of flux and reveals the property of flux that is being 
measured when the line interface tension is measured. 

It will be shown that the effect of the boundary conditions is to reweight flux con- 
figurations according to the iV-ality of the flux through an interface. With z\ = e 2m ^ N , 
configurations with iV-ality n get a factor of z n = z\. Thus the partition function with 
interface boundary conditions is the sum over n of sets of configurations with iV-ality n and 
with an additional factor of z n in the weight. This is unlikely to be relevant to a physical 
situation. I do not know how this reweighting would arise in a natural way. 

When matter is added, the situation changes significantly. Matter fields in the funda- 
mental representation explicitly break the global Z(N) symmetry. The N line states are 
no longer degenerate. Above the critical temperature, there could be metastable states of 
the line variables separated by walls jnj. However, it follows from the matter-free analysis 



that there are no physical metastable states or walls in the presence of matter. A physical 
effect of the matter is to lower the free energy. This comes about through an increase in the 
entropy when new configurations with sources are allowed. 

In this paper, I will use a very simple Z(2) model to discuss these issues from first 
principles and in detail. The more interesting SU(N) and Z(N) gauge theories can be 
handled with the same approach, but they require more theoretical elaboration. I will cover 
those in another paper. 

While calculations in terms of the physical flux variables and in terms of the unphysical 
lines lead to the same result, the associated interpretations are completely different. Inter- 
faces in the line variables may be a convenient device for making approximate calculations 
of physical quantities in terms of unphysical variables. However, one should be cautious 
with heuristic arguments that rely upon the physical reality of the interfaces. In terms of 
physical variables, there are no interfaces or metastable states. 

Section II covers gauge fields without matter. It is shown that there are no interfaces 
in the physical variables associated with the interfaces in the lines. Boundary conditions 
designed to detect the contribution from line interfaces are applied in cylinder geometry. 



When the problem is expressed in the physical variables, it becomes clear that there are no 
associated physical interfaces. The associated physical effect is the "conservation" mod N 
of the total TV-ality of the flux along the cylinder axis, which is perpendicular to the line 
interface. The Z(2) model is used to illustrate these points. Section III introduces matter 
of nontrivial iV-ality. First, there is a short argument that because there are no physically 
distinct high-temperature phases of the Yang-Mills field, there are also no metastable states 
in the presence of matter. Matter simply lowers the free energy. It is shown that this can 
be understood in terms of the physical variables as an increase in entropy. Section IV is a 
brief conclusion. 

II. NO PHYSICAL WALLS IN HOT YANG-MILLS THEORY 

This section treats gauge theory without matter fields. The Hamiltonian formulation 
in the A = gauge is used. This is the most physical formulation since it is closest 
to the elementary language of Hamiltonian quantum mechanics. A lattice cutoff is used. 
The physical degrees of freedom are gauge group elements on links of the spatial lattice. 
One basis of states specifies a definite value of the group element on each link. The other 
convenient basis specifies that the amplitude to be at different points in the group is given 
by a matrix element from an irreducible representation. This is flux on links. Either is 
a description of physical variables. Physical states are configurations of physical variables 
with the additional condition that the Gauss law constraints are satisfied. To enforce these 
constraints, the unphysical variables are introduced in the Lagrangian formulation. These 
are the group elements on links in the fourth direction or the fourth component of the gauge 
field. 

If the partition function is written as a sum over physical states, then the unphysical 
variables have no work to do and need not even appear. Clearly, interfaces in the unphysical 
variables are not physical per se. So the question is: do they reflect the existence of interfaces 
in the physical variables? 
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One possibility is that interfaces in the Wilson lines are associated with gaps in the 
percolating flux. There could be large islands of flux that have no flux links to the rest of 
the lattice. However, such an island can be connected to the rest of the flux network with 
the local addition of a finite number of links. For large islands, there are many ways to do 
this. Thus, the probability of a large island appearing is very small — much smaller than 
the probability for a wall to appear in the line variables. This follows from the fact that a 
local, finite-energy change in the line variables does not destroy the wall, it only moves its 
position a little. This conclusion can also be deduced from the calculations below. 

Previous work on the high-temperature, homogeneous phases of the lines |3|-|5| makes 
the existence of physical interfaces seem very unlikely. It was shown that the physical 
configurations associated with the N pure phases of the Wilson lines at high temperature are 
the same. Thus, in terms of physical variables, there is a single high-temperature ensemble. 
There are not N physical phases that could be separated by walls. 

So there does not appear to be any reason why there should be interfaces in the physical 
variables. In particular, they are not required by symmetry as in the case of the unphysical 
variables. The walls in the line variables do not imply the existence of walls in the physical 
variables. The following discussion reveals the property of physical flux that is associated 
with the unphysical interfaces in the lines. It is the "conservation" mod N of flux along the 
spatial direction perpendicular to the interface. 

The following discussion deals mainly with the Z(2) flux model. This model is equivalent 
to the Z(2) spin Ising model 0. It is also an approximation to Z{2) gauge theory. Thus, 
the discussion and conclusions also apply to Z(2) gauge theory. 

The gauge and flux models are related through the basis states and because of the 
limited role of the Hamiltonian at high temperature. In the Hamiltonian formulation of the 
3+1-dimensional gauge theory, a basis for states is flux on links. More precisely, on each 
link, an irreducible representation matrix element for the group is specified. The allowed 
configurations are restricted to those that satisfy the Gauss law constraints. In the case of 
Z(2), the flux is either zero or one on each link. The constraints are that each site must 
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have an even number of links with ones. These are exactly the configurations of the flux 
model. 

In the flux model, 9(1) = 1 if there is flux on link I, and it is if there is not. The energy 
of a link of flux is a. Configurations are weighted with a Boltzmann factor 

e~ H/T with H = J2<j9(1). (2.1) 

i 

In a gauge theory, the Hamiltonian comes from the transfer matrix. There is a factor of this 



form [[Ty] and a factor from the spatial plaquette term in the Lagrangian. The plaquette 
factor is diagonal in the basis where the links have definite group elements, but it is not 
diagonal in the flux basis. Thus it could complicate the discussion that follows. However, it 
is sufficient to consider the high temperature limit. This is the region where the Wilson lines 
are ordered. It is the large (3 limit of the spin system. It is the region where the interfaces 
in the lines are sharp. In a theory, such as the Z(2) case, with discrete, finite states, the 
Hamiltonian is a perturbation at high temperature. Thus, the flux weighting of ( [2.1| ) plays 
no significant role in the discussion. That does not change if a plaquette term is added to 
H. Thus, the discussion in this paper is sufficient for high-temperature, Z(2), gauge theory. 

It is important that there is a discrete and finite set of flux states for each link. It might 
appear that a new approach will be needed if there is an infinite number of irreducible 
representations as there is for any continuous group. However, it turns out that the role 
that is played by flux in the Z(2) case will be played by iV-ality for the SU (N) gauge groups. 
Quite a bit more machinery is required, but the basic structure of the development is not 
changed. 

The analysis begins with the Z(2) Ising flux model used in previous work ||. In this 
model, space is a three-dimensional cubic lattice. There are variables 9 on links that can 
have the values 1 or to indicate the presence or absence of flux. These are the physical 
variables. A configuration is specified by the function 9(1). States with definite flux \9) are 
labeled by that function. A general state has a wave functional that gives the amplitude for 
the system to be found in the various basis states of definite flux. The energy of a link with 



flux is a. The weight for a configuration is 

{0\e- H ' T \6) = e -^' a0{l) = He-T* 9 ®. (2.2) 

i 

The sum over configurations is restricted to those in which the number of links at a site that 
have flux is even. Let the collection of all such configurations be C . It is a subset of the 
unrestricted collection of all configurations C . The partition function is 

Z = Tr'\e- H ' T } = *E(9\e- H ' T \9) = $>4 (2.3) 

c c 

The equivalence of this flux model to the Ising model results from using site variables to 
enforce the restriction on configurations, and then doing the 9 sums. The site variables are 
the Ising spins. First, consider the sum of 9(1) over the 2d = 6 Vs contained in the set J(z) 
of links with endpoint v. 

m = e m (2-4) 

To force this to be even at each site, introduce the site variables s(i) that take the values 
±1. The factor 

\ E *« S(i) (2-5) 

s(i)=±l 

has the desired effect. A factor like this is introduced into the partition function sum for 
each site i. The spins s(i) are the same as the Wilson lines of the Z(2) gauge theory. This 
model is equivalent to the Ising model for the spins s(i) with the Ising (3 related to the 
gauge T by e~ a l T = tanh/3. Note that I never refer to the Ising model 1/(3 as temperature. 
Temperature always refers to the T appearing in (|2.2|) . Large T gives large (3, and ordered 
Ising spins. It is convenient to introduce /i = a/T. In the high T region, 

e~ w ~ fi/2 (2.6) 

To discuss interfaces, it helps to introduce boundary conditions that force at least one 
interface in the lines to appear. Use solid-cylindrical geometry: space is finite in the trans- 
verse x and y directions and very long in the z direction. There are periodic boundary 
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conditions in the x and y directions. The transverse area in dimensionless lattice units is 
A. Let L be large and positive. Apply the boundary conditions s = — 1 for z = —L/2 and 
s = 1 for z = L/2. The limit of interest is L — > oo followed by A — > oo. 

First review the picture in terms of the unphysical spin variables. The high-temperature 
region is the large-/?, ordered phase for the Ising spins. The interface is approximately flat 
and has a Boltzmann weight e~ 2/3A . Thus, the energy per unit area a is related to f3 by 
a/T PS 2(3 as (3 — > oo. With L — 1 planes between z = —L/2 and z = L/2, there are L 
places to put an interface. The boundary conditions require an odd number of them. Since 
j3 is large, a is large, and the interfaces are sparse. The sum in the partition function can 
be extended to infinity to obtain 

2Z> = e Le ~ aA/T - e~ Le - aA/T . (2.7) 

The partition function Z with the same boundary conditions at each end is 

2Z = e Le - aA/T + e- Le ~ aA/T . (2.8) 

The ratio is 

Z'/Z = e^ F '- F ^ T = 1 — (F' — F)/T — 1 — 2e' 2Le "" A . (2.9) 
The excess free energy 

F' — F = 2Te~ 2Le ~ 2/3A (2.10) 

decays to 1 exponentially at a rate determined by the interface energy per unit area a = 2f3T. 
This relates the activity for the wall e~ 2l3A to Z'/Z and F' — F. The approximations are 
valid for large L and Le~ 2/3A 3> 1 . 

Now consider the same situation in terms of flux variables. In the high-temperature, 
deconfined phase with /i = a/T small, there is dense, percolating flux. The flux is almost 
random except for the constraint that there be an even number of links with flux at each 
site. The following discussion shows how the boundary conditions, which require at least 
one interface, affect the flux configurations contributing to the partition function. 
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The system consists of all the links and sites strictly between the two planes with the 
fixed spins. To enforce the constraints, there is a factor of 

Z s 

at each site. For a site on the z = —L/2 boundary, this becomes (— \) e . The longitudinal 
links that end on boundary sites are not constrained on the boundary end because the spin 
is fixed there. Of course, the constraints at the sites on the other ends of those links must 
still be satisfied. The boundary links can be freely chosen to satisfy the constraints on the 
first layer of spins at each end. On the positive z end, nothing else happens. However, on 
the negative z end, there is a factor of —1 for each link with flux coming into the system 
from the boundary. 

From the local Gauss law constraints, it follows that the total flux Y^, X y^z{x,y,z) mod 
2 on a transverse slice of longitudinal links is independent of z. This will be referred to as 
flux "conservation". Thus each configuration can be described as even or odd. The effect 
of the boundary conditions is to weight the odd configurations with and extra factor of —1. 
The partition function sum 

Z= J2e~ H / T + J2^ H/T (2.12) 

even odd 



is replaced by 



Z'= Y,e~ H / T -Y,e~ H/T - (2-13) 

even odd 



This gives 

Z'/Z = (Z e - Z )/(Z e + Z ) = 1- 2e-^~ F ^ T . (2.14) 

The last step is correct for large (F a — F e ). I will show in a moment that this free energy 
difference is proportional to the length L, so that ( |2.14| ) is consistent with ( |2.9| ). Thus, the 



free energy difference per unit length of odd flux verses even flux is the quantity of interest. 

The next part of the discussion shows the direct connection between the spin interfaces 
and the conservation of flux. Consider the structure of the partition function. The role of 
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the spins is to enforce the local constraints on flux. In the calculation leading to ( |2.9|) , the 
spins are constant on transverse planes. Thus, they are enforcing the weaker constraint of 
flux conservation. The details follow. 

Let the constant value of the spins at position i along the longitudinal axis be S(i). In 
this situation, the flux on transverse links is unconstrained. In the partition function, the 
transverse flux sum gives the same factor on each plane. That factor is raised to the power 
L — 1. This is an uninteresting, overall factor. Another ingredient in the partition function 
is the flux sum on a single slice of longitudinal links. There are A longitudinal links into a 
transverse plane. The sum becomes a product of A independent sums and gives 

£ e -^ [S(i)S(i + = [1 + e-»S(i)S(i + 1)} A = a + bS(t)S(t + 1) (2.15) 

6 

with 

a=h(l + e-<*) A + (l-e-») A } and b = h(l + e^) A - (1 - e^) A \. (2.16) 

Now consider the calculation of Z' . Let G(i) = J2j ^ e ^ ne sum °f a ^ ^ ne longitudinal 

flux on a slice at position i. 

Z ' = E ^(-l) e ( 1 ) e ^ ( 1 )»S(l) e(1)+0(2) e-^ e(2) 5(2) e ( 2 ) +e ( 3 ) . . . (2.17) 
e s 

e - M e(7V) s ^ N s0(N)+e(N+i) e -tie(N+i) _ (2.18) 

First do the 9 sums to obtain 

Z' = - bS(l)][a + bS{l)S{2)] . . . [a + bS{N - l)S(N)][a + bS{N)}. (2.19) 

s 

This brings us back to a one-dimensional Ising model. The spin sums can be done. The 
results are 

Z' = a L - b L (2.20) 

and 

Z'/Z = 1 - 2e- Lln(a/b) . (2.21) 
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For small /i, a/b ~ 1 + 2(^) A . Comparing ( |2.21| ) with ( |2.14j) gives 



Z'/Z = 1 - 2e~ 2L( t) A (2.22) 

and 

(F -F e )/T = 2L(^) A . (2.23) 



Since /i and (3 are related by ( |2.6|) this is equivalent to ( |2.9| ). 

This shows that the interfaces in the unphysical spin variables are associated with flux 
conservation i.e. the total flux mod 2 on a transverse slice of links is independent of longi- 
tudinal position. 

With that understood, one can redo the calculation of Z' by a much easier method using 
flux variables. The partition function Z' is the difference of two contributions. It is the sum 
of all the even configurations minus the sum of all the odd configurations. The weight for 
one slice of even longitudinal flux is 



£i5>-/*V» = a. (2.24) 



For odd longitudinal flux, it is 



^ = E^E^ E V^ +1 = &- (2.25) 

For small fi, the configurations of flux on longitudinal links are almost independent from 
slice to slice except that the total flux mod 2 is the same. Thus 

Z' = Z e -Z D = z L e -z L a =a L -b L (2.26) 

as in (|2.20|) . This is almost zero because the weights for even and odd configurations are 
almost the same when /i is small. 

This calculation of Z' is simple and direct in the physical flux variables. There is no 
reference to interfaces. The result is the same as that from the spin calculation, but the 
associated physical picture is completely different. In flux, it is the Gauss law constraints 
rather than interfaces that are important. 
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This completes the main discussion of the matter-free case. Before considering the in- 
troduction of matter in the next section, let us examine the approximations that have been 
used. The main approximation has been to assert that after summing over flux configu- 
rations on transverse links, the configurations on longitudinal links from slice to slice are 
independent except for the restriction that they have the same total flux mod 2. This is 
equivalent to the statement that at small fx and large f3, the spins on a transverse slice are 
ordered. 

To examine the approximation, fix a configuration of flux on longitudinal links subject 
only to the constraint of flux conservation. Consider a plane of transverse links separating 
two slices of longitudinal links. In the large T, fx — > limit, the allowed configurations get 
equal weight. The density of flux is high. Without the constraints, it would be 0.5 per link. 
The constraints say that only configurations with an even number of links with flux at each 
site are allowed. If there were no longitudinal flux, then that constraint would apply to the 
transverse flux on the plane. The effect of the longitudinal flux is to require odd transverse 
flux at some sites. Since longitudinal flux is conserved, there are an even number of such 
sites on each slice. The crucial fact is that the number of allowed configurations of transverse 
flux is independent of the number of pairs of odd sites. Given that, it is clear that the sum 
over all transverse flux configurations is independent of the longitudinal configuration. This 
is the approximation that was used above. 

The problem is reduced to showing that the number of allowed configurations is inde- 
pendent of the number of pairs of odd sites. Consider two configurations of pairs of odd 
sites c and c'. The corresponding sets of allowed flux are C and C . I need to show that 
the numbers of configurations in C and C are the same. In spin variables, this is just the 
statement that n-point functions of spins are position independent at infinite (3. In the flux 
picture, the result can be established by producing a one-to-one mapping between C and 
C. Without loss of generality, assume that c and d differ by c' having a single extra pair of 
odd sites relative to c. Connect those two sites with some path of links. Reversing the flux 
<-> 1 on the links of the path gives the desired map. 
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The conclusion of this section is that with the cylinder geometry and boundary conditions 
that reveal the presence of walls in spin variables, there are no walls in physical flux variables. 
The spin walls are associated with flux conservation along the cylinder axis. The ratio Z' jZ 
of partition functions with odd and even boundary conditions can be used to extract the 
surface energy of a wall. In the flux picture, it is a measure the difference in free energy per 
unit length between configurations with even flux and those with odd flux on the longitudinal 
links. 

However, the universe does not have this artificial geometry and boundary conditions. 
For this free energy difference to be physically relevant, it would require that a factor of — 1 
be associated with a closed surface when the flux through it is odd. This was contrived in 
the discussion above through the choice of geometry and boundary conditions. I do not see 
how it could arise naturally. Thus, the spin interface tension has a physical counterpart but 
it is not physically relevant. 

III. EFFECTS OF MATTER 

With the introduction of matter carrying an irreducible representation of SU(N) or 
Z(N) with nontrivial A-ality, the global Z(N) symmetry is explicitly broken. This lifts the 
degeneracy of the multiple, high-T phases of the line variables. When the effect from the 
matter is small, there will be metastable states. Such states and their possible astrophys- 
ical consequences have been discussed [|Kj. The important effect of the matter is that it 
breaks the global symmetry. This can be studied in a simpler situation where the matter is 
represented by an external field h linearly coupled to the Wilson lines: h*L + hL*. 

In a general situation with Z(2) global symmetry, the constraint effective potential with 
h 7^ is a double well with one side lower than the other. As h — > 0, this becomes the 
symmetric double well with degeneracy. Thus the existence of the metastable state is linked 
to the existence of degeneracy. This is the generic situation away from the critical point in 
the two-phase region. 
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However, without matter, there is a single, physical, high-temperature phase. Thus, 
there is no degeneracy that could be connected to multiple states, some metastable, when 
symmetry breaking is present. I conclude that in the physical flux picture, degeneracy, 
interfaces, and metastable states are all absent. 

Nevertheless, matter does affect the free energy. It is interesting to describe this effect 
entirely in the flux picture. The idea is to see what is happening to the flux configurations 
when one side of the effective potential for the lines is lowered. Lowering the free energy 
means increasing the partition function. I will show that this increase in Z comes from 
an increase in the number of allowed flux configurations in the presence of matter i.e. an 
entropy increase. 

It is sufficient for this discussion to consider the case of heavy, dilute matter. It follows 
from the work in Ref. top |12| that this is well approximated by a small external field. Thus, 
in the Z(2) flux model, the symmetry breaking effect of matter will be represented by a 
weak external field h coupled linearly to the line or spin variables. 

In the high-temperature, two-phase region, the free energy is lowered. When the partition 
function is expressed in terms of flux variables, the effect of h is to allow odd sites with a 
density controlled by h. I will show that this increases the number of allowed configurations, 
increases Z, and decreases F. 

The extra factor that comes from the introduction of matter (as represented by the 
external field) is 

e E t hs(i) = -Q e hs(i) = -Q( cosh h + g ^ sinh ft) (3 J) 

i i 

The sums on the spin variables give the constraints on flux. For each site, there is a factor 
of 

Se cosh h + 5o sinh h. (3.2) 

The shorthand 8(e,o) stands for the restriction that the flux at the site be (even, odd). With 
h small, most sites have even flux, but there is a small density of sites with odd flux. The 
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sites with even flux have a slightly reduced weight cosh/i vs. 1, while sites with odd flux are 
now allowed with a small weight sinh h. I will show that the latter effect is more important 
so that F decreases while Z increases. 

A little notation is convenient. Let Z(ii, . . . , i n ) be the partition function from summing 
over configurations with odd flux at the sites i\, . . . , i n and even flux at all others. In finite 
volume V, n must be even. Further, define Z(ii, . . . , i n ) so that it vanishes if n is odd or if 
any two i's are equal. Let Z(h) be the partition function with the external field h. 

Z(h) =Y,e~^^ ° m J] [cosh hd E + sinh h6 ] (3.3) 

i 

= Z(0)(cosh/i) y [l + ^(tanh/ i )"^ £ Z(h, . . . , i n )/Z(0) (3.4) 

n n - ii,...,i n 

Since only even n's are present, Z(h) is even in h —>■ —h. 

It is useful to extend the definition of Z to odd n. Factorization properties of Z are also 
needed. It will be shown that there is a function Z(ii, . . . ,i n ), defined for even and odd 
n with the following properties. For even n, Z(i±, . . . ,i n ) = Z(i\, . . . ,i n ). To describe the 
factorization, let P = {pi,P2, . . .} be a partition of {ii, . . . , i n }. When the clusters of points 
are widely separated, 

z(i u ...,i n )^Ylz{ie Pl ). (3.5) 

i 

For example, 

Z{ii,i2) = Z(i 1 ,i 2 ) — > Z(ii)Z(i 2 ) as \i x - i 2 \ — >• oo (3.6) 

and 

Z(h, i2i is) —> Z(ii, i2)Z(is) for |z 3 — ii| — > oo with i\ and i 2 fixed. (3.7) 

Note that for odd n, Z is not a finite-volume partition function with an odd number of 



odd-flux sites. Z(h) is defined by replacing Z with Z in (|3.4j ). The original Z(h) can be 
recovered via 

Z{h)= l -\Z{h) + Z{-h)}. (3.8) 
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It is a bit more convenient to deal with 



Z' = l + ]T(tanh/iri- £ Z{i 1 ,...,i n )/Z(0). (3.9) 



n! „ 
ii,...,i n 



Now introduce by 



Z(h) = e w{h) (3.10) 



and write 



OO -| 

iy(/i) = ^(tanh/i) ri - X! W{h,...,i n ). (3.11) 



, n! . 

n=l ij, ...fa 



W(i\, . . . ,i n ) is defined recursively from Z(ix, . . . , i n ) as usual |T3",n|. If Z(ii, . . . , i n ) 
factorizes for large separation, then W(ii, . . . , i n ) goes to zero if any separation between z's 
is large. Thus, the sum on the i's in ( p.ll|) is proportional to the spatial volume V. The 



quantity of interest is the coefficient of V in the h — >• limit. This comes from w = W(i). 
Combining Q, (H), (|3~l0| ), and (|3~TT|) gives 



= ^(O)^^)^-^ 1 ^ 1 + e - vlhwl ] (3.12) 

for V large and /i small. In this limit, the e y][hw]{ term dominates. Thus as h increases from 
zero, there is an increase in Z and a decrease in the free energy. It is clear from the discussion 
that this is a direct consequence of the increase in entropy that comes from including some 
configurations with odd-flux sites. This effect dominates the middle factor (cosh/i) y , which 
reflects the change in the weight of the configurations with even flux at every site. This 
completes our discussion of the behavior of the system in terms of the physical variables. 
However, the conclusions are based on properties of Z and Z that must now be established. 



A. Properties of Z{i\, . . . , i n ) 

It must be shown that Z has factorization properties which allow Z to be defined and to 
factorize according to (|3~5|). First recall that T is large and the flux is dense. I will use an 
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expansion in y(l) = e~^ a6 ^ — 1, which is small when T is large. Let C(i\, . . . , i n ) be the set 
of flux configurations with odd flux at the indicated sites and even flux at all others. The 
number of configurations in C is denoted by j^C. Clearly, n must be even. 

Z(h,...,i n )= £ Y[e~^ (3.13) 

C(ii,...,i n ) I 

= J2Ui 1 + y( 1 )] (3-i4) 

= Et 1 + E?/( / i)+ E y(h)y(k) + ■ ■ ■] (3.15) 

h {h,h} 

= E[i + EA E y('i)---y(W] (3-i6) 

m m - l 1 ...,l m \l i ^U 



An m-fold / sum in ( [3.15| ) is over all subsets of the links with m elements. In Q3.16|) , the 



/-sums are restricted only by the condition that no two links can be the same. Introduce 
the notation z(j) and w(j) for the generating functions of the ?/(/)'s. 

*U) = e u ® = ^I$ + 'L^ E y(h)...y(i m )j(ii)...j(l m )} (3.17) 

it C m m - l u ...,l m 

^[i+E^T E <h,...,lm)j(h)...j(lm)] (3.18) 
m m - h,...,l m 

and 

OO 1 

»0') = E-f E w(h,...,i m )j(h)...j(i m ). (3.i9) 

m=l " b - h,...,l m 



The dependence of z's and w's on the z's has not been shown here. Equation Q3.18| ) is a 



definition of the z(l\, . . . , l m ). The k sums are unrestricted in ( |3.18| ) so z(l\, . . . , Z m ) vanishes 
if any Vs coincide. 

The key technical result is a factorization property of z(li, . . . , l m ). Let R be a partition 
of the links . . . ,l m } into subsets r,. Whenever the links of each r, are well separated 
from those of all other rj, then 

«({iei2}) = n^er i }). (3.20) 

i 

This will be proven below. 
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Given fl3.20|) , the usual combinatoric arguments lead to the conclusion that w(h, . . . ,l n ) 
falls to zero rapidly if any link separation becomes large. This means that the sums in ( |3.19| ) 
can be considered as a sum on all possible clusters of links followed by a sum that moves the 
cluster to all possible positions. Since every link that is included brings a factor of y ~ 
only small clusters contribute at high temperature. 

Consider the case in which the odd sites are partitioned into groups pi that are widely 
separated. The objects of interest are 

W U = 1 ) = Y.— 7 w(iep 1 ,iep 2 ,...;l 1 ,...,l m ) (3.21) 

ill I — 

m ln - h,...,l m 

and 

ln Z(yH ! 7f ' n ^ n ^ = E~j E ep 1 ,iep 2 ,...;l 1 ,...,l m )- w(l x , l m )\ (3.22) 
Z(uj m ml h ^ lm 

= w(pi) + w(p2) + . . . + w{resi). (3.23) 

w(pi) is the contribution from the part of the I sums where the cluster of links {Zi, . . . , l m } 
is around the group of sites Pi. w(rest) is from the terms in the I sums where the cluster 
of links is far from any odd site. When the cluster of links {Zi, . . . , l m } is far from an odd 
site i, then the flux configurations around the cluster, which determine the value of w via 
( p.l7| )- (|3.19| ), are independent of the position i. This is because C(i,j, . . . , k) can be mapped 
to C(i',j, . . . , k) without changing the configurations away from i and near j, . . . ,k. This 
can also be shown using the technique introduced at the end of this section. Thus, the last 
term w(rest) is independent of i. 

So, Z(ii, . . . , i n ) is a product of factors — one for each cluster pj of odd sites. The number 
of sites in a cluster can be even or odd. These factors can be used to define Z(ix, ■ ■ ■ ,i n ) for 
odd n. So defined, the Z(i±, . . . , i n ) have the required properties. 

Equation ( |3.23| ) shows that 



W{i) = w(i). (3.24) 

It is straightforward to calculate w(i) in low orders of the high-temperature expansion. The 
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first nonzero term comes at order T 6 when the site i has contributing y's on all the links 
attached to it. The result is 

e w& = 1 _ _L (i_ e -"/*y (3.25) 
32 

= 1 - 2e~ l2f3 (3.26) 

= (s(i)) (to this leading order). (3.27) 



The the right hand sides of ( |3.26| ) and ( p.27 ) could also be obtained from a simple calculation 



in terms of spin variables. This shows that the flux calculation produces the same result as 
the spin calculation. 

Finally, the factorization of z in its I variables shown in (|3.20|) must be derived. This 
requires several steps and a little notation. The configurations have odd flux at sites . . . , i n 
only. This is understood and not explicitly displayed in the following discussion. Subsets of 
links are indicated by X, X\, etc. A specific configuration of #(Z)'s, I G X is denoted by G. 
The set C(Oi, 6 2 , . . .) is all allowed configurations of 6*'s with given, fixed sub configurations 
on the specified subsets of links and with odd sites at ii, . . . ,i n . 

The derivation establishes the following assertions: 

1. There is a one-to-one map C(O) <-» C(O'). 

2. Thus, #C(6) = #C{&). 

3. If two subsets of links X\ and X 2 are separated, then a configuration of #'s on X 2 is 
allowed or not independent of the configuration of the #'s on X\. 

4. #c(e 1 ,e 2 ) = #c(e 1 ,6' 2 ) 

5. There is a factorization 

(/i(e 1 )/ 2 (e 2 )) = (/ 1 (e 1 )>o(/ a (e 2 )> . (3.28) 



The subscript indicates an ensemble in which all the contributing configurations are equally 
weighted — the T = oo limit. The relation ( |3.20|) is a special case of (|3.28| ). 
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Here is the reasoning that establishes each of the assertions. 

1. Consider a configuration Co in which the subconfiguration on X is 6 and another 
configuration c' with 6' on X. For c , let F be the set of links on which 9(1) = 1, and 
for c' , the corresponding set is Y' . The map is defined as a reversal (0 <-> 1) of the #'s 
on links in U = (Y U Y') — Y n F'. This is an invertible map on and onto C(ii, . . . , i n ) 
that sends c — > c' . To see that the same map, with c and c' unchanged, sends any 
configuration in C(0) to one in C(O'), consider links I E X. If 0(7) = 1 in and ©', 
then / e Y fl 7' and I G" U, and 6>(/) does not change. If 9(1) = in 6 and 6', then 
/ £ F U F' and / £ U, and 0(Z) does not change. If 9(1) = 1 in and 0(Z) = in 6', 
then I G F and / ^ F' and / G Z, and 0(/) is reversed, as it should. Clearly, the fourth 
case works like this one. 

2. This is an immediate consequence of the fact that the map is one-to-one and onto. 

3. Let X be a cube of links which contains X 2 . Let R be a larger cube of links containing 
X. The region outside R, which contains X 1 is V , and the region in R but not in 
X is /. For any configuration on V Ul, I can change links in I so that there is 
either no path of 9 = 1 links or one path of 9 = 1 links that connects X and Vq. This 
depends on whether the number of odd sites in X is even or odd. Now change the 
configuration on X 2 to any other allowed configuration. With appropriate changes in 
I and no changes in Vo, it can be connected to Vq. Thus, an allowed configuration on 
X 2 is allowed independent of the configuration of 0's well separated from it and, in 
particular, those in X±. 

4. With X 1 and X 2 well separated, take X = X x U X 2 . The previous steps combine to 
show that there is a one-to-one map C(©i,© 2 ) C(Qi,Q 2 ). Thus, #C(©i,© 2 ) = 

#C(©!,© 2 ). 

5. Consider 
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(A(©i)/ 2 (© 2 ))o = ^77 E /i(©i)/a(©2) (3-29) 

= ^EE E A(©i)/ 2 (© 2 ) (3.30) 

" ei e 2 c(e 1 ,e 2 ) 

^E/i(9i)E/ 2 (e 2 ) E (3.3i) 

ei e 2 c(ei,e 2 ) 

= {^E1E][ E ]}(/i(©i)>o(/ 2 (© 2 ))o. (3.32) 
" ei ©a c(8i,e 2 ) 



The first factor {. . .} in ( |3.32|) is one, and the relation is established. 



In this section, I have shown that the effect of fundamental representation matter can 
be understood in physical variables. The usual picture in terms of unphysical variables is 
that the Z(2) symmetry is broken, and the double well of the effective potential becomes 
asymmetric with one side lowered in free energy and the other raised to a metastable state. 
In terms of physical variables, there is no double structure with a metastable state. The 
free energy is simply lowered. This results from an increase in the entropy that comes from 
allowing some new configurations where the gauge field does not transform as a singlet at 
the sites with matter. 



IV. CONCLUSION 

In terms of the Wilson line variables, the high temperature phase of a gauge theory is 
associated with broken symmetry. This is a multiphase region. Phase separation and domain 
walls are possible. However, in a Hamiltonian formulation, the lines are unphysical variables 
so it is interesting to formulate the results in terms of physical variables and physical states. 
I have done this for the Z(2) case. In terms of the physical variables, there is only one high 
temperature phase. There are no interfaces. The interface tension of the lines has a physical 
counterpart. It arises when configurations with odd flux are given an extra factor of —1 in 
the sum for the partition function. It is hard to see what physical situation might make this 
modification relevant. 



22 



When matter that breaks the center symmetry is added, there is a lowering of the free 
energy of one of the high temperature phases of the lines. The other phases are then 
metastable. In the physical variables, there are no metastable states. However, the free 
energy of the single high-temperature phase is lowered by the matter. The origin of that 
is an increase in the entropy that follows from allowing some new configurations where the 
gauge field does not transform as a singlet at the sites with matter. 
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